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Abstract 

ABSTRACT 

In this paper we study the critical behavior of the two-dimensional antiferromagnetic Ising model 
in both uniform longitudinal (H) and transverse ($7) magnetic fields. Using the effective-field theory 
(EFT) with correlation in single site cluster we calculate the phase diagrams in the H — T and 
£1 — T planes for the square (z = 4) lattices. We have only found second order phase transitions 
for all values of fields and reentrant behavior was not observed. 
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I. INTRODUCTION 



In the last decade there has been an increasing number of works dealing with magnetic 
models to study quantum phase transitions. In particular, considerable interest has been 
directed to the transverse Ising model (TIM) used to describe a variety of physical systems []]- 
1 It was originally introduced by de Gennes fl pseudospin model for hydrogen-bonded 
ferroelectrics such as KH 2 PO^ in the order-disorder phenomenon with tunnelling effects so, 
it has been successfully used to study a number of problems of phase transitions associated 
with order-disorder phenomena in other systems {4], [5]. It provides a good description to 
analyze some real anisotropic magnetic materials in a transverse field. 

Theoretically, various methods has been employed to study the criticality of the TIM 
such as renormalizaton group (RG) method [6], effective field theory (EFT) [t q| , mean 
field theory (MFA) [lOj, cluster variation method (CVM) [11], pair approximation (PA) 12), 
Monte Carlo (MC) simulations [13] . and so on. The critical behavior of the one dimensional 
TIM has already been established through exact results, where the ground-state energy, 
the elementary excitations and the correlation functions were obtained [l^]. The TIM is 
among the simplest conceivable classes of quantum models in statistical mechanics to study 
quantum phase transition [l5, [lq . The ferromagnetic TIM has been studied intensively. 
The critical and thermal properties of the ferromagnetic and antiferromagnetic TIM are 
equivalent, but the properties of these two model are very different at not null longitudinal 
field (H 7^ 0). For example, the ground-state phase transition in the ferromagnetic TIM is 
smeared out by the longitudinal field in contrast to the antiferromagnetic TIM for which 
the phase diagram remains qualitatively the some at (H ^ 0). 

Experimental investigations on the metamagnetic compounds such as FeBr2 and FeCl 2 



17] and Ni{NO^) 2 2H 2 18|, 119] . under hydrostatic pressure, have been performed. For 



example, in the FeBr 2 compound was observed a sharp peak in magnetization measurements 
under a field inclined by 33° with the c axis (perpendicular to the plane) of the crystal. They 
concluded that the peak was affected by the ordering of the planar spin components. It is 
obvious that the field can be decomposed into the longitudinal and transverse fields. The 
model is described by the Ising Hamiltonian to which is added a term which represents 
the effects of the transverse field part. Due to the requisite of non-commutativity of the 
operators in the Hamiltonian, deriving the eigenvalues of the Hamiltonian is a very difficult 
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problem. Therefore, many theoretical methods have been used to investigate this system 
20 

From a theoretical point of view it is known that the effect of the transverse field in 
the Ising model TIM is to destroy the long-range order of the system. Many approximate 
methods |20h29| have been used to study the critical properties of this quantum 

model. Some years ago, a simple and versatile scheme, denoted by differential operator 
technique [30|, was proposed and has been applied exhaustively to study a large variety of 
problems. In particular, this technique was used to treat the criticality of the TIM 3l| 
obtaining satisfactory quantitative results in comparison with more sophisticated methods 
(for example, MC). This method is used in conjugation with a decoupling procedure which 
ignores all high-order spin correlations (EFT). The EFT included correlations through the 
use of the van der Waerden identity and provided results which are much superior that 
the MFA. The TIM was first studied by using EFT 3l| for the case of spin S = 1/2, and 
generalized [25] for arbitrary spin-S* ^1/2. 

On the other hand, the critical and thermal properties of the transverse Ising antiferro- 
magnet in the presence of a longitudinal magnetic field have been few studied in the literature 



33j. Using the c 



assical approach MFA and the density-matrix renormalization-group 
method (DMRG) [33j, the ground state phase diagram in the (H — Q) plane was studied 
in one-dimensional lattice. Critical line separates the antiferromagnetic (AF) phase with 
long-range order (LRO) from the paramagnetic (P) phase with uniform magnetization. The 
quantum critical point S c = (j) c decreases as h c = (-j) increase, and is null at h c = 2.0. 
The MFA approach does not give the correct qualitative description of the critical line 33]. 
Firstly, the quantum fluctuations shift the ground critical point S c = 1.0 to S c = 2.0 at 
h = underestimates critical value. Second, the form of the critical line shows an incorrect 
behavior around of the critical point h c = 2.0. 



By using the EFT approach, Neto and de Sousa [34( have studied the ground-state phase 
diagram of this quantum model on two-dimensional (honeycomb (z = 3) and square (z = 4)) 
lattices, where was discussed the possibility of existence of a reentrant behavior around 
h c = z critical value. The spin correlation effects are partially taken into account in EFT, 
while it is entirely neglected in MFA. The differences of results given by using EFT and 
MFA show that the spin correlation has important on the phase diagram. The ground-state 
phase diagram in the (H — Q) plane is qualitatively similar to the results of Fig. ([T]) in Ref. 
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33[, but the reentrant behavior found by MFA occurs near h c = z. 

In the present paper, using EFT we investigate the quantum phase transitions of the Ising 
antiferromagnetic in both external longitudinal and transverse fields. This work is organized 
as follows: In Sec. II we outline the formalism and its application to the transverse Ising 
antiferromagnetic in the presence of a longitudinal magnetic field; in Sec. Ill we discuss the 
results; and finally, in Sec. IV we present our conclusions. 



II. MODEL AND FORMALISM 



The model studied in this work is the nearest-neighbor (nn) Ising antiferromagnet in a 
mixed transverse and longitudinal magnetic fields divided into two equivalent interpenetrat- 
ing sublattices A and B, that is described by the following Hamiltonian 

h = j E - H E < - n E °f » c 1 ) 

<i,j> i i 

where J is the AF exchange coupling, denoted the sum over all pairs of nearest-neighbor 
spins (z) on a ci-dimensional lattice (here we treat the two-dimensional lattices with z = 3 
and 4), erf is the v[= x,z) component of the spin-1/2 Pauli operator at site i, and H(Q) is 
the longitudinal(transverse) magnetic field. 

The erf and erf spin-1/2 Pauli operators do not commute, then a nonzero field (Q) trans- 
verse to the spin direction causing quantum tunneling between the spin-up and spin-down 
eigenstates of erf and quantum spin fluctuations. These fluctuations decrease the critical tem- 
perature T c at which the spins develop long-range order. At a critical field Q c , T c vanishes, 
and a quantum phase transition between the AF ordered state and a quantum paramagnetic 
state occur. To the best of our knowledge, the model ([T]) at finite temperature (T 7^ 0) has 
not yet been examined in the literature. In particular, the ground-state phase diagram was 
studied by Neto and de Sousa [34|, then, in this work we generalize it to analyze the field 
effects at finite temperature by using EFT. 

The ground-state of the model (CD) is characterized by an antiparallel spin orientation 
in the horizontal and vertical directions and so it exhibits Neel order within the initial 
sublattices A and B (see Fig. ([1])), that is denoted by the AF state. 

The competition between the antiferromagnetic exchange field presents interesting prop- 
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FIG. 1: Ground state of the quantum Ising antiferromagnet on an square lattice described by 
Hamiltonian given in Eq. (fT|). 



erties in the phase diagram. In particular, the model (IIJ) has an AF (ordered) phase in 
the presence of a field, with the decreasing transition temperature as the fields intensity 
increases, where at T = (ground- state) a second-order transition occurs at both critical 
field with H c (= zJ) and Q c . 

To treat the model (0Q) on two-dimensional lattices by the EFT approach, we consider 
a simple example of cluster on a lattice consisting of a central spin and z perimeter spins 
being the nearest-neighbors of the central one. The nearest-neighbor spins are substituted 
by an effective field produced by the other spins, which can be determined by the condition 
that the thermal average of the central spin is equal to that of its nearest-neighbor ones. 
The Hamiltonian for this cluster is given by 



where A and B denote the sublattice. 

From the Hamiltonians fl2]) and fl3]), we obtain the average magnetizations in sublattice 
A, rriA = (cia)> an d B, m# = (crfg), using the approximate Callen-Suzuki relation {3]], that 
are given by 





and 




(3) 
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'"^ = < ^— . 4 , „, tanh/V(i/-„ ,l 2 <J- ) . ID 

and 

m B = { —^~"\ B tanh/V(ff-a 1B ) 2 + f] 2 ), (5) 



( ^ = tanh /3 y/(H-a 1B )* + w) , 

\ y/(H - ais) 2 + / 

z 



z z 

where a 1A = JJ2 a (i+s)B and a iB = JJ2 a (i+s)A- 

Now, using the identity exp(aD x )F(x) = F(x + a) (where D x = is the differential 
operator) and the van der Waerden identity for the two-state spin system (i.e., exp(acrf) = 
cosh(a) + erf sinh(a)) the Eqs. (BJ and ([5]) are rewritten as 



m A 

and 



X\{a x + a z {1+5)B P x )) F{x)\ x= ^ (6) 



with 



™ B = ( Y[(a x + <t z {1+s) M ) F(x)\ x=0 , (7) 

Wo / 

Fix) = , H ~ % = t&nh[3J(H -x) 2 + n 2 , (8) 

where a x = cosh(JD x ) and /3 X = sinh(JD x ). The Eqs. (jSJ) and ([7]) are expressed in terms 
of multiple spin correlation functions. The problem becomes unmanageable when we try to 
treat exactly all boundary spin-spin correlation function present in Eqs. fl6]) and (J7|). In this 
work we use a decoupling right-hand sides in Eqs. ([6]) and (EJ), namely 

(<riA<7jB ■ ■ ■ a iA) - m A m B ...m A , (9) 
where i ^ j ^ . . . ^ I and = (of) (// = A, B). The approximation fl9]) neglects 
correlation between different spins but takes relations such as ^(of^) 2 ^ = 1 exactly into 
account, while in the usual MFA all the self- and multispin correlations are neglected. Using 
the approximation ((9]), the Eqs. (jSJ) and (J7j) are rewritten by 

z 

m A = Y,MTN,H,Q)m p B , (10) 

p=0 

and 

z 

m B = J2MTN,H,Q)m p A , (11) 

p=0 



with 



A P (T N , H, fl) = —^—a>-*% F(x)\ x=0 , (12) 

where the coefficients A p (Tn,H,Q) are obtained by using the relation exp(aD x )F(x) = 
F(x + a). 

In terms of the uniform m = Jmi+mj) and staggered m s = |(m,A— m#) magnetizations, 
and near the critical point we have m s — > and m — > mo, the sublattice magnetizaton tua 
expanded up to linear order in m s (order parameter) is given by 

m A = X (T N , H, tt, m ) + Xi(T N , H, tt, m )m s , (13) 

with 

z 

X (T N , H, Q, m ) = M T n, H, Q)m p , (14) 

p=0 

and 

2 

X 1 (T N , H, n, m ) = - Y,pMTn, H, fi)m^ 1 . (15) 

p=0 

On the other hand, in this work only second-order transitions are observed, therefore, to 
study the phase diagram only we analyze the Eqs. (fl4"|) and (fT5!) in the limit of m s — > 
on can locate the second-order line and using the fact that tua = rriQ + m s in Eq. ( flBl we 
obtain 

X Q (T N , h, 5, m ) = m , (16) 

X 1 (T N ,h,6,m ) = l, (17) 

at the critical point in which m s = 0, 5 = Q/J and h = H/J. 

Thus, we can get an analytic solution for second-order transition where m s is the order 
parameter which is used to describe the phase transition of the model ([T]). The magnetiza- 
tions of two sublattices are not equal for m s ^ 0, and the system is in the antiferromagnetic 
phase. The magnetizations of two sublattices are equal for m s = 0, and the system is in the 
saturated paramagnetic phase. 
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III. RESULTS AND DISCUSSION 



The numerical determination of the phase boundary (second-order phase transition) is 
obtained by solving simultaneously Eqs. ( I16p and (fT7|) . We determined the phase diagrams 
of this quantum model in the h — T and 5 — T planes for the square (z = 4) lattice that 
comprises a field-induced AF phase (m s 7^ 0) at low fields and a P phase (m s = 0) at 
high fields. In the case h = 0, we have = —ms {m = ) a nd the critical behavior 
reduces to the transverse Ising ferromagnetic analyzed in Ref. 3l|. For z = 3 and z = 4, 
the values of the quantum critical point for h — obtained are 5 C = 1.83 (<5^ FA = 3.00) 
and 8 C = 2.75 (5f FA = 4.00), respectively |34|. In the limit of null fields h = 5 = 0, we 
have obtained the value m = and kBT N /J = 3.085 (square lattice) and ksT^jJ = 2.104 
(honeycomb lattice). The results for the critical temperature obtained by simple EFT with 
cluster of N = 1 spin (EFT-1) ksT^/J = 3.085, when compared with the exact solution 
^bTn/J = 2.269 . . ., for the Ising model on a square lattice, is not quantitatively satisfactory. 
On the other hand, with the increases of the size cluster (N = 2, 4, 9, . . .) we have a possible 
convergence of the results for the critical temperature. In the present paper we have only 
interested in obtaining qualitative results for the phase diagram. 

The study of quantum phase transition is nowadays one of the main areas of researching 
in condensed matter physics, so most the experimental and theoretical studies are developed 
to magnetic quantum critical point (QCP). For the square lattice (z = 4), the QCP for h = 
obtained by using EFT [3lJ, 5 C = 2.75, can be compared with the results found by using other 
methods, for example, S c = 4.00 of MFA, S c = 3.22 of path-integral Monte Carlo simulation 
36], S c = 3.05 of the DMRG |37| . S c = 3.08 of the high-temperature series expansion 1|5], 



5 C = 3.02 of the effective-field renormalization group (EFRG) |38j, S c = 3.00 of the PA |l2|, 
and 5 C = 3.04 of the cluster MC [sqJ. With the increases of the size cluster we will have a 
possible convergence for the results rigorous of 5 C 37-391. 

In Fig. (J2J), the phase diagram in the h — T plane is presented for the square lattice, 
with selected values of 5. In this figure, the order of the phase transition between the AF 
and P phases is invariably of second-order for all values of the transverse field. The critical 
temperature decreases monotonically as the longitudinal field (H) increases, and so is null 
at h — h c (5). the critical behavior of h c versus transverse field (5) (ground-state phase 
diagram) has been analyzed by Neto and de Sousa [34 1. 
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The same qualitative phase diagram is also observed in the 5 — T plane in Fig. [3] for 
the square lattice, with selected values of h. When the h longitudinal field increases the 
critical curve T c (8) decrease. We show that there is no reentrant behavior in the region of 
low temperature, and also the phase transition is of second order for all values of the fields 
h and S. 

We hope that the present method can be employed for more complex models, for instance, 
the generalization of the model to treat three-dimensional lattice, where multicritical points 
in the phase diagrams are observed in MFT 4o|, 41]. This reentrant behavior does not occur 



at low temperature and that was found in a classical approach (MFA) that does not give 
the correct description of the ground-state phase diagram, with the presence of reentrant 
behavior near h = h c and the quantum critical point 5 c /z = 1.00 is independent of the 



coordination number z 



33 



34]. 



IV. CONCLUSIONS 



We have used the effective-field theory with correlation in single site cluster to obtain the 
state equations of the TIM antiferromagnetic. We study the phase diagrams in h — T and 
5 — T planes for the square (z = 4) lattice. The transverse field 5 has important effects on the 
phase diagrams because it destroys the long-range order of the system. The results show that 
for a given 5, the critical longitudinal magnetic field h decreases with increasing temperature. 
In the phase diagrams there are no reentrant phenomena. Our results (EFT) are consistent 
with second-order transitions from the antiferromagnetic phase to paramagnetic phase at 
null transverse field. The qualitative results for the phase diagrams can be obtained by using 
more reliable methods such as quantum Monte Carlo simulation and renormalization group 
approaches. Furthermore, the investigations of this three-dimensional model are expected 
to show many characteristic phenomena. They will be discussed in the work. 
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FIG. 2: Dependence of the reduced critical temperature ksT/J as a function of the reduced 
magnetic field h = H/J for a square lattice and several values of 5. The results are qualitatively 
identical to the honeycomb lattice. 
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FIG. 3: Dependence of the reduced critical temperature ksT/J as a function of the reduced 
magnetic field 6 = Q/J for a square lattice and several values of h. The results are qualitatively 
identical to the honeycomb lattice. 
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